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ABSTRACT 

We present Atacama Large Millimeter Array (ALMA) Cycle-0 observations of the CO J = 6 — 5 
line in the advanced galaxy merger Arp 220. This line traces warm molecular gas, which dominates 
the total CO luminosity. The CO emission from the two nuclei is well resolved by the 0.39" x 0.22" 
beam and the exceptional sensitivity and spatial/spectral resolution reveal new complex features in 
the morphology and kinematics of the warm gas. The line profiles are asymmetric between the red 
and blue sides of the nuclear disks and the peak of the line emission is offset from the peak of the 
continuum emission in both nuclei by about 100 pc in the same direction. CO self-absorption is 
detected at the centers of both nuclei but it is much deeper in the eastern nucleus. We also clearly 
detect strong, highly redshifted CO absorption located near the southwest side of each nucleus. For 
the eastern nucleus, we reproduce the major line profile features with a simple kinematic model of a 
highly turbulent, rotating disk with a substantial line center optical depth and a large gradient in the 
excitation temperature. The red/blue asymmetries and line-to-continuum offset are likely produced 
by absorption of the blue (SW) sides of the two nuclei by blue-shifted, foreground molecular gas; 
the mass of the absorber is comparable to the nuclear warm gas mass (^10^ Mq). We measure an 
unusually high Tco/Afir ratio in the eastern nucleus, suggesting there is an additional energy source, 
such as mechanical energy from shocks, present in this nucleus. 

Subject headings: galaxies: ISM, galaxies: active, techniques: interferometric, galaxies: kinematics 
and dynamics 


1. INTRODUCTION 

As the closest example {D = 77 Mpc, 1"= 373 pc) of 
an ultra-luminous infrared galaxy (ULIRG), Arp 220 has 
been studied intensively at a wide range of wavelengths. 
It is an advanced merger of two gas-rich galaxies with 
high dust optical depth (rioo/xm ^ 5; Rangwala et al. 
2011) and a large reservoir of molecular gas 10^^ 
Mc^) concentrated in its nuclear region (jScoville et al.l 
I1997D . It hosts one of the most extreme star-forming 
environments in the local universe and is a widely used 
template for interpreting observations of high- 2 :, dusty 
star-forming galaxies. High-resolution ground-based ra¬ 
dio interferometric observations sho w that it has two nu¬ 
clei separated by rou ghly 1" (e.g.. I^oville et all 119971: 
iSakamoto et al.l l2QQ9f ) , with evidence of a highly dust- 
obscured AGN in its w estern nucleus (jPownes fc EckarO 
I2QQ7I: [Engel et al.ll2Qllf) . Additionally, a kiloparsec-scale 
central molecular disk surrounding the two nuclear disks 
has been inferr ed from CO J = 2 — 1 interferometric 
observations by iScoville et al.l (|1997D . 

Observations of high-J CO rotational transitions 
from the Herschel Space Observatory of nearby galaxies 
(including Arp 220) and Galactic star-forming regions 
have determined that there are generically two com- 
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ponents in the molecular gas: a warm (high-pressure) 
component traced by mid-to-high-J lines that comprises 
a minority of the mass (^10% of the total molecular gas 
mass) but dominates the CO luminosity (^90% of the 
total luminosity), and a cold (low-pressure) component 
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PapadoDoulos et ^12014 ). Comparison of the observed 


CO spectral line energy distributions (SLEDs) and 
far-infrared luminosities to theoretical models indicates 
that in several prominent cases (e.g., Arp 220, M 82, 
NGC 6240, M83) the gas heating is dominated by some 
mechanism other than UV or X-ray photons (as in 
photon-dominated and X-ray dominated regions) or 
cosmic rays; mechanical energy from shocks generated 
by supernovae and stellar winds are the most plausible 
sources. A recent study by Greve et al. (2014) using 
CO measurements in 62 local (U)LIRGs and 35 high-z 
galaxies also concludes that mechanical energy is the 
more likely source of energy for this warm component. 
The pressure in the warm component is found to be 
typically two orders of magnitude larger than the cool 
component pressure (e.g., 10^ K cm“^ vs. 10^ K cm“^). 
These findings suggest that this warm component of 
the molecular gas is intimately associated with the 
feedback processes, which play an essential role in 
galaxy formation and evolution. 

Our ongoing Herschel archival survey is systematically 
deriving physical conditions of the molecular gas in a 
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large sample (^200) of nearby galaxies by modeling the 
combined observations from the Herschel SPIRE-FTS 
instrument, which provides CO SLEDs from J = 4 — 3 
to J = 13 — 12, and ground-based measurements of low- 
J CO lines. The results from the analysis of a n ini¬ 
tial sample of 17 galaxies (jKamenetzkv et aDl2Q14l ) con¬ 
firmed the main findings described above. This study 
also found that the mid-J CO transitions (J > 4 — 3), 
where the CO spectral line energy distribution peaks, 
primarily arise in the warm component and their line 
luminosities are well correlated with the total CO lumi¬ 
nosity. Thus these transitions can be used as a reliable 
tracer of warm molecular emission. Because of its rela¬ 
tively large beams, Herschel could not spatially resolve 
most galaxies, so only average properties of the warm gas 
have been inferred, and it has been difficult to associate 
the warm gas with the distribution of star formation, 
cool molecular gas, and nuclear activity. 

The CO J = 6 — 5 line is now accessible with ALMA, 
which can finely resolve the morphology and kinematics 
of the warm gas in many nearby galaxies. This transition 
has been previously obs erved in Arp 220 with the SMA 
(|Matsnshita et al']l2QQ9[) . However, the lesser resolution 
and lower sensitivity of SMA at these higher frequencies 
could not resolve the emission from the two nuclei. In this 
paper, we present the first high-resolution and high-S/N 
maps of the warm molecular gas in Arp 220 using the CO 
J = 6 — 5 line. The corresponding continuum observa- 
tions (at 435/im) were presented in IWilson et al.l (|20l4 
Paper I hereafter). These observations showed that the 
western nucleus is considerably warmer and more opti¬ 
cally thick (Tdust = 200 K,r 434 ^rn = 5.3) than the east¬ 
ern nucleus (Tdust = 80 K,r 434 ^rn = 1-7). The two nuclei 
combined account for roughly 80% of the total infrared 
luminosity of Arp 220. 

The ALMA observations and data reduction are pre¬ 
sented in Section 2. The observed morphology and kine¬ 
matics in the CO J = 6 — 5 emission is described in 
Section 3 followed by detailed modeling of the kinemat¬ 
ics of the nuclear disks in Section 4. We discuss disk 
stability and the energetics of the warm molecular gas in 
Sections 5 and 6, respectively. Detection of absorption 
in the SiO J= 16 — 15 is presented in Section 7 and our 
conclusions are summarized in Section 8. 

2. OBSERVATIONS AND DATA REDUCTION 

Arp 220 was observed with ALMA on 2012 December 
31 (Project 2011.0.00403.S) at the end of cycle-0 in an ex¬ 
tended/hybrid configuration. The Band 9 receivers were 
tuned to cover the CO J=6-5 transition (rest frequency 
691.473 GHz; Arp 220’s redshift = 0.01813) and adjacent 
continuum, with a total bandwidth of 8 GHz; however, 
the usable bandwidth was somewhat reduced due to over¬ 
lap between the spectral windows covering the GO line 
and the need to flag a handful of end channels in each 
of the four spectral windows. The observations covered 
projected baselines from 13 to 374 m. The correlator 
was configured in low-resolution wide-bandwidth TDM 
mode to give a spectral resolution of 15.625 MHz. The 
total on-source integration time was 44.4 minutes using 
23 antennas. 

The data were flagged, calibrated, and imaged by 
ALMA staff prior to delivery. The data have two unusual 
technical problems. First, no absolute flux calibrator was 


observed as part of the data set. The absolute flux scale 
was set using a flux of 6.0 Jy for the bandpass calibrator, 
3G279, determined by a Band 9 observation taken within 
one day of these observations. This process introduces 
some additional uncertainty on the absolute flux scale 
which is difficult to quantify. Second, the signal-to-noise 
ratio on the phase calibrators is very low, which required 
a major adjustment to the standard cycle-0 calibration 
process. The process used by ALMA staff was as fol¬ 
lows. (1) Obtain a bandpass calibration solution from 
the bandpass calibrator. (2) Use the bandpass calibrator 
to solve for phase offsets among the four spectral windows 
and two polarizations. (3) Obtain a solution for phase 
versus time from the bandpass and two gain calibrators 
while combining the data for both polarizations and all 
four spectral windows together. (4) Obtain a solution for 
amplitude versus time from the three calibrators again 
combining the data from all polarizations and spectral 
windows. Use this solution to estimate the fluxes (0.32, 
0.40 Jy) of the two secondary calibrators. (5) Obtain a fi¬ 
nal solution for amplitude versus time using the adopted 
flux of the two gain calibrators. 

For our analysis, we began with the delivered cali¬ 
brated uv data and made new, cleaned data cubes fol¬ 
lowing the same steps carried out by ALMA staff using 
GAS A version 4.1.0. Repeating the imaging allowed us 
to adjust the clean boxes to reduce the level of artifacts 
in the final cubes even further from the very good images 
originally delivered to us and to improve the resolution 
by switching to uniform weighting. We also imaged all 
four spectral windows and adjusted the velocity ranges 
used to make the continuum images (jWilson et al.ll2014 ) 
to avoid all emission and absorption features. 

The reduction notes from the ALMA staff noted that 
self-calibration was required to get the best images. Fol¬ 
lowing their method, we carried out two rounds of phase- 
only self-calibration using natural weighting and a rela¬ 
tively small clean box to include emission from the bright 
central region. A final round of self-calibration on both 
amplitude and phase significantly improved the residuals 
in the data cube. As noted by ALMA staff, we see signif¬ 
icant amplitude corrections are required and these cor¬ 
rections correlate among the different spectral windows. 
These steps significantly improved the noise level in the 
cube, from 36 mJy/beam to 16 mJy/beam in the line 
free channels and from 100 mJy/beam to 18 mJy/beam 
in channels where the line emission was strong. 

Before making the final image cube, we carried out 
continuum subtraction using a zeroth order fit to most 
of the channels in the two spectral windows located away 
from the bright GO emission. In the final imaging step, 
the clean box was adjusted with frequency to include 
all regions with significant emission, which improved the 
negative residuals in certain channels with bright emis¬ 
sion. We made image cubes with both natural and uni¬ 
form weighting to maximize sensitivity and angular res¬ 
olution, respectively. We used a velocity resolution of 
10 km s~^ for the GO data cube. As a final step, the 
uniform cube was smoothed slightly so that all channels 
would have the same beam (0.39" x 0.22") at position 
angle 28°. The typical la noise in the cleaned uniform 
cube is 30 mJy beam“^ or 0.9 K per channel. 

After a small correction for primary beam attenua¬ 
tion, the total integrated GO J = 6 — 5 flux is 3670 
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Figure 1. Channel maps covering the CO emission and absorption over ~800 km s The 1% contour level for the continuum emission 
is shown relative to the peak continuum flux of 1.38 Jy beam“^. Redshifted CO absorption can be seen in the high velocity channels (5505 
- 5705 km s“^). Note: The maps and spectra in this paper are not corrected for the small primary beam attenuation. 


Jy km This value is in agreement with the to¬ 

tal flux obtained from Hersche l-SFIRE of 4070 ± 80 Jy 
km s“^ (jRangwala et al.ll2011l ) within their calibration 
uncertainties (about 10% for Herschel-SPIRE and 15% 
for ALMA). Even in the absence of an absolute flux cali¬ 
brator, the final flux calibration of our data cube is con¬ 
sistent with literature values. 

The channel maps from the uniform data cube are pre¬ 
sented in Eigure[T]for CO emission over ^800 km s“^. 
The extent of the continuum emission is indicated by the 
1% contour level. Deep CO self-absorption is detected 
around the systemic velocity of the eastern nucleus that 
appears to go below zero. Also, redshifted CO absorp¬ 
tion is seen in higher velocity channels. These features 
are seen in the spectra at the continuum positions even 
before the continuum is subtracted. There are gaps in 
the data cube over two velocity intervals: 4245 - 4305 
km s“^ and 5735 - 6035 km s“^ (see Eigure El). The 
latter gap is close to the redshifted CO absorption and 
limits our ability to measure the full velocity extent of 
this absorption feature. These features are further dis¬ 
cussed in B3. 


3. OBSERVED MORPHOLOGY AND KINEMATICS OF THE 
WARM MOLECULAR GAS 

The CO J = 6 — 5 emission traces the warm molecular 
gas peaks at the centers of the eastern and western nuclei, 
which are unambiguously resolved in the 0.39" x 0.22" 
resolution CO map (Eigure [2^), and has an extended 
component surrounding the two nuclei. This overall mor¬ 
phology a£rees_with__thnCOJ=J— and CO J = 3 —2 
maps of iSakamoto et all (j 19991 . l2008[) , in which the emis¬ 
sion associated with the two nuclei was characterized as 
compact nuclear disks that are embedded in an extended 
outer disk. A simple Gaussian fitting to the two nuclei 
shows that the western nucleus is almost circular with 
a deconvolved size (EWHM) of ^ 0.55" (205 pc using 
1"= 373 pc), while the eastern nucleus is elongated in 
the NE-SW direction with a deconvolved size (EWHM) 
of 0.74" X 0.59" (276 pc x 220 pc). This indicates that 
the CO emission is extended compared to the dust con¬ 
tinuum emission, for which the half-maximum diameters 
are 76 x (< 70) pc and 123 x 79 pc for the western and 
eastern nucleus, respectively (Paper I). However, there is 
likely extended, low surface brightness dust emission that 
is below the current sensitivity of our observations: the 
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Figure 2. ALMA high-resolution images of the CO J = 6 — 5 line and 691 GHz continuum in Arp 220. (a) CO integrated intensity 
(moment-0) map made from uniform cube. The 10%, 20%, 40%, 60% and 80% contour levels are shown relative to the peak flux of 389 
Jy beam“^ km s“^ in the western nucleus, (b) Mean velocity Field (moment-1), (c) Velocity dispersion (moment-2) map. The 20%, 40%, 
60% and 80% continuum contours are shown relative to the peak continuum flux of 1.15 Jy beam“^ in the western nucleus, (d) Peak 
temperature (moment-8) map and integrated CO intensity contours in White, (e) The peak of the line emission (color map) is offset from 
the peak of the continuum emission (solid contours) by roughly lOOpc. The SiO J = 16 — 15 line absorption (dashed contours - 20%, 40%, 
60% and 80% of -70 Jy beam“^ km s“^) is detected only in the western nucleus and is unresolved, (f) Moment-0 contours of CO J = 6 — 5 
redshifted absorption over 5505 - 5725 km s“^ (yellow contours). The continuum emission (color) and line emission (solid contours) are 
overplotted. The yellow contour levels are 20%, 30%, 50%, 70% and 90% relative to -25 Jy beam“^ km s“^. The contours depict both the 
compact and extended components of the redshifted absorption. 


combined continuum emission at 692 GHz detected in the 
two nuclei is about 80% of the total continuum emission 
measured in the ALMA map (Paper I) and only 56% of 
the t otal continuum emissio n measured in the Herschel 
data (jRangwala et alJl2Qllf ). The peak flux, integrated 
flux and source sizes for CO J = 6 — 5 emission are listed 
in Table 1. 

The overall CO J = 6 — 5 velocity field (Figure [JJ)) 
is also similar to the CO J = 2 — 1 and CO J = 3 — 2 
maps, which show a large velocity gradient in the NE - 
SW direction. The velocity extent is ^500 km s“^ within 
0.3" of the nu clei, again similar to the CO J = 2 — 1 
velocity extent (|Sakamoto et al.iri999[ ). The CO J = 6—5 
velocity dispersion (cr = FWHM/2.3548) map (Figure 
[21(c)) shows the dispersion peaks are coincident with the 
continuum peaks, with values as large as 175 km s“^ 
and 160 km s~^ for the eastern and western nucleus, 
respectively. 

To comp are with the CO J = 2 — 1 and CO J = 3 — 2 
maps from [Sakamoto et al.l (|2008l ). which have a beam 
FWHM of 0.5", we convolved the CO J = 6 — 5 line map 


to match their beam size. Using these maps we find that 
about 50% of the total detected CO J = 6 — 5 emission 
arises from within the nuclear region (i.e., within a radius 
of 0.45"). This fraction is higher than for the low-J lines 
(^30%), implying that the warm molecular emission is 
relatively more compact. 

The peak temperature (moment-8) map. Figure [21(d), 
displays two noticeable features in comparison with the 
integrated flux map: the semicircular ridge of emission in 
the western nucleus is more prominent in peak temper¬ 
ature than in integrated flux, and there is an off-nucleus 
“hot spot” in the eastern nucleus (although its tempera¬ 
ture is lower than the nuclear peak). Comparable high- 
resolution maps in additional CO transitions will be re¬ 
quired to determine whether these temperature features 
correspond to variations in excitation temperature. 

Overplotting the continuum contours on the CO J = 
6 — 5 integrated intensity map (Figure [2^) shows an offset 
between the peak of the continuum emission and the peak 
of the line emission in both nuclei; this offset is about 0.2" 
and has the same direction. This offset is not seen in the 
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Table 1 

ALMA Arp 220 line Measurements 


Quantity 

Eastern Nucleus 

Western Nucleus 

Units 




GO J = 6 — 5 Emission 

Peak flux location 

Integrated flux (1" x 1") 

Peak flux 

Gaussian FWHM 

Gaussian-flt Integrated flux 
Line-to-Continuum ratio 
Lco/Lpm^ 

Deconvolved FWHM 

Warm Gas Mass^ 

Warm Gas Golumn Density‘s 

(15:34:57.30, +23:30:11.53) 
1070 

330 

0.83 X 0.63 

1565 

1035 

6.71^4 X 10-^ 

0.74 X 0.59 

1.4 X 10® 

9.2 X 10^2 

(15:34:57.22, +23:30:11.68) 
1305 

515 

0.66 X 0.63 

1795 

990 

Q.otfe X 10-® 

0.59 X 0.54 

1.6 X 10® 

1.5 X 10^3 

J2000 

Jy km s“^ 

Jy beam“^ km s“^ 

// 

Jy km s“^ 
km s“^ 

// 

M© 

H/cm-2 

1 SiO J = 16-15 Absorption 

Peak flux location 

Peak flux 

Gaussian FWHM, PA 
Gaussian-Fit Integrated flux 
Mean depth 

Velocity extent 

No Detection 

(15:34:57.22, +23:30:11.43) 
-80 

0.385 X 0.255, 27 

-87 

-0.23 

>400 

J2000 

Jy beam“^ km s“^ 
(",°) 

Jy km s ^ 

Jy beam“^ 
km s“^ 


Note. — The peak and integrated fluxes corrected for primary beam attenuation are reported from the natural 
weighted (higher sensitivity) moment-0 map with a beam size of 0.447" x 0.273" (PA = 28°), while the source 
FWHM is reported from the higher-resolution uniform weighted map with a beam size of 0.39" x 0.22". 

The S/N in our data is very high - ~55 in the peak flux and > 400 on the integrated flux - leading to very 
small measurement errors; thus they are not reported. The overall uncertainties are dominated by systematic 
calibration uncertainties of ^^15%. 

^ Using the LpiR measured from the ALMA 691 GHz continuum emission in Paper 1. Note that the large range 
in uncertainty in this ratio arises entirely from the large systematic uncertainty on the LpiR. 

^ Using M(iL2)warm/-Lc'06-5 = 21 Mq/Lq derived from Herschel FTS CO SLED (Rangwala et al. 2011) 

° Using the Gaussian FWHM and Area = tt 
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Figure 3. Continuum subtracted spectral line profiles at the lo¬ 
cation of the continuum peak for the eastern and western nuclei, 
in the reference frame of CO. 

low-J lines (up to J = 3 — 2) from previous observations. 
We explain in Section 4 that this offset is only apparent, 
and is likely caused by absorption of the entire SW side 
of the line emission in the two nuclei. 

Moreover, CO absorption is detected at the centers of 
the two nuclei as shown in Figure (H the absorption is 
much deeper in the eastern nucleus and goes below zero 
in the continuum subtracted spectra. This absorption 
is also shown in the channel map close to the systemic 
velocity of the eastern nucleus (Figure [U channel 5345 


km s“^). The central dip seen in the eastern nucleus sug¬ 
gests the presence of absorption; that this dip goes below 
zero (these are continuum subtracted spectra) indicates 
both that this feature is the result of absorption and that 
we are seeing absorption of the 6 — 5 line against the con¬ 
tinuum source, i.e., the absorbing gas is absorbing both 
line and continuum emission. 

The deep CO absorption at the center was not reported 
from the interferometric observations of the CO J = 2 —1 
and 3 — 2 lines with 0.5" resolution ([Sakamoto et al.l 
20081) . In Figure [H the CO J = 3 — 2 profile from 
Sakamoto et al.l (| 2 0091 . private communication 2015), ob¬ 
served using a 0.3" beam, is overplotted on our CO 
J = 6 — 5 profile (at the continuum peak). This com¬ 
parison clearly shows that for the eastern nucleus a sim¬ 
ilar absorption dip is present in the low-J CO l i ne, al - 
though it was not discussed by iSakamoto et al.l (|2009f ) 
(further discussion is provided in section 4). In fact, the 
CO J = 3 — 2 line shows a double-peaked profile in the 
eastern nucleus similar to that of CO J = 6 — 5, except 
the blue side is much weaker than the red side. Compar¬ 
ison of the western nucleus profiles also shows similarity 
between the two transitions, i.e., a central absorption dip 
and an overall complexity in the line profiles. 

The spectrum of the western nucleus in Figure [3] (bot¬ 
tom panel) also shows redshifted CO absorption at ve¬ 
locities between roughly 5505 km s~^ and 5725 km s“^, 
much higher than the systemic velocity of the system. 
Absorption is also detected towards the eastern nucleus 
but in off-center positions and in higher velocity channels 
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(see Figure [T] (channel 5705 km s“^) and Figure [5|). The 
redshifted absorption is shown with yellow contours in 
Figure [2jf), which also shows the CO emission in solid 
contours and the continuum emission in color. The ab¬ 
sorption has a deep, compact component as well as a 
shallow, extended component. We note that although 
the shallow extended component is present over many 
velocity channels, the less than ideal calibration of the 
data and the low S/N mean that we cannot be certain 
it is not an artifact of the cleaning process. The peak 
of the redshifted CO absorption is also offset in the two 
nuclei relative to the continuum peak but in the oppo¬ 
site direction compared to the line emission peak, i.e., 
towards the SW direction. Because this absorption is 
present towards both nuclei, its plausible that it arises 
from a continuous feature covering the southern part of 
the whole system, i.e., at least 400 pc in length (roughly 
the distance between the two nuclei). This continuous 
redshifted absorption feature could be interpreted as an 
infalling molecular filament (see section 4.2), not unex¬ 
pected in a chaotic late-stage merger. It is difficult to 
ascertain the full velocity extent of this feature because 
of gaps in our data but it is at least 200 km s~^ wide. 

A highly excited SiO line, J = 16 — 15, is detected 
in absorption (dashed contour in Figure Efe); also see 
Figure [3]), but only towards the western nucleus. This 
line is discussed further in Section 7. 

In summary, the exceptional sensitivity and spectral 
resolution from ALMA has allowed us to detect several 
components in the warm molecular gas in Arp 220 show¬ 
ing more complexity in the morphology and nuclear kine¬ 
matics than previously observed. 



Velocity [radio, LSR](km/s) 


Figure 4. Compa rison between the CO J = 6 — 5 (this work) 
and CO J = 3 — 2 (jSakamoto et al.l 120091 . private communication 
2015) velocity profiles at the location of the continuum peak for 
the eastern and western nuclei in 0.39" x 0.22" and 0.3" beams, 
respectively. 


4. MODELING: KINEMATICS OF THE NUCLEAR DISKS 


Both the eastern and western nuclei show clear kine¬ 
matic evidence for rotating disks. In addition to the ve¬ 
locity gradients clearly seen in Figure [Sfb), in Figure [5l 
we show spectra along the major axes of both nuclear 
disks, spaced at multiples of a half-beamwidth (^80 pc) 
from the dust continuum peaks. In both cases we see 
clear signs of rotation: a systematic shift in the velocity 
of the emission peak with position along the major axis, 
with an isolated red or blue peak on opposite sides of the 
continuum peak. The kinematics of the eastern nucleus, 
in particular, appear relatively simple, with a double- 
peaked (but asymmetric) profile at the continuum cen¬ 
ter that shifts to single red- or blue-shifted peaks with 
increasing distance along the major axis. In comparison, 
the velocity profiles of the western nucleus are more com¬ 
plex. For this reason we present here a simple toy model 
intended to explain the main features of the spectra for 
the eastern nuclear disk, and then discuss the implica¬ 
tions for both nuclei and the larger-scale gas distribution 
around the nuclei. 

In the presence of velocity dispersion and systematic 
motion (e.g., rotation or outflow), we can write the line 
optical depth as a function of frequency as 

Tv = To(f>v ( 1 ) 

where tq is the line center optical depth and the Doppler 
line profile (for non-relativistic velocities) is given by 


(pj, = exp 


Aud 


{uo/c)n-v 

Aud 


( 2 ) 


where Aud is the Doppler width (v^ times the dis¬ 
persion), z^o is the line center frequency, n is the di- 
rection of propagation, and v is the velocity (e.g., 
lHumme~ fc Rvbickilll968[ ). 

The nuclear disk is substantially inclined to our line of 
sight, with i 50° based on the deconvolved continuum 
source size (Paper I). For simplicity, we assume that the 
disk is perfectly edge-on. We also assume that the disk is 
uniform in the vertical direction, which means that the 
disk thickness is only a meaningless scaling parameter. 
This reduces the model to two dimensions, which con¬ 
siderably speeds up the numerical modeling. We expect 
that the inferred model parameters would only change 
by factors of order unity in going from this edge-on, 2D 
disk to a tilted, 3D disk, e.g., the true disk rotation curve 
will increase by 1/cosi from the model curve. The only 
features not captured by this approximation would be 
systematic motions that are confined largely or entirely 
to the disk vertical direction. While our simple edge-on 
model can only include features that are a function of 
projected distance from the disk center with no depen¬ 
dence on angle within or height above the disk, the match 
between the model spectra and the data are sufficiently 
good (see 4.2) that features violating this restriction are 
not required to explain the observed spectra. 

We choose a coordinate system such that the x-axis lies 
in the plane of the sky, the ^-axis is now the disk vertical 
direction, and the ^-axis lies in the disk and is oriented 
along the line of sight towards the observer, with 0 the 
angle in the disk with respect to the x-axis (so 6> = 7r/2 
points towards the observer). The radial coordinate r 
then lies in the plane of the disk, and the distinction be- 
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tween cylindrical and spherical radius has been removed 
(e.g., a spherical outflow lies purely in the disk plane in 
this approximation). The projected line-of-sight velocity 
is then 

n-Y = Vc{r) cos 9 + Vout{r) sin 9 (3) 

for circular rotation velocity Vc and purely radial outflow 
velocity Vout- The distortion of the line profile function 
by the systematic velocity terms plays a key role in deter¬ 
mining the appearance of the disk emission. Note that 
equation m with equation ([3j) is antisymmetric about 
^ = 7r/2 in the rotation-only case, and symmetric in the 
outflow-only case (provided that the outflow is itself sym¬ 
metric about this line), which cuts the needed amount of 
computation in half. However, it is neither in the case 
that both rotation and outflow are present. 

We specify the line and continuum source functions 
and optical depths as a function of radius in the disk, 
which is assumed to be circular and axisymmetric. Both 
the source functions and the optical depth profiles are 
assumed to be Gaussians in radius rI3 With both line and 
continuum emission and absorption, the source function 
Sjy is given by the ratio of the combined emissivities 
to the sum of the absorption coefficients 


Su 


ju,l T ju,d 


( 4 ) 


with the I and d subscripts referring to line and dust, 
respectively. The optical depth at frequency u is now 


Tu = Tu,l + Ti,^d- (5) 

In terms of the line-only and dust-only source functions, 
equation (j4]) can be written 


S. 


T GLv^dSv^d 
T ^u,d 


(6) 


i.e., the individual source functions weighted by their re¬ 
spective absorption coefficients, or, equivalently, by the 
corresponding values of . Equations and dni) can be 
simplified by noting that we can take the dust emissiv- 
ity and absorption coefficient (and thus the dust source 
function) to be independent of frequency over the width 
of the line profile. 

The emitted intensity along a line of sight through the 
disk is then 


lu 




a^dz . 


( 7 ) 


To solve this equation, we divide the disk into segments 
along the z—direction at each value of impact parame¬ 
ter p (perpendicular distance from disk center); changing 
variables from distance along the line of sight to optical 
depth riy{z) = aj^{z')dz'^ is then the sum over 

pieces of the form 


Su{z){l - (8) 


^ Note that this is not a physical model, in the sense that the 
excitation of the CO line has not been calculated from conditions 
in the disk: the line source function and optical depth are specified 
independently of one another, and both the line and continuum 
source functions have arbitrary normalizations. This model is in¬ 
tended to aid in understanding the main features of the emission 
from the eastern nucleus, not reproduce it in detail. 


where Tiy{z) is the local optical depth, and tI^^{z) is the 
total optical depth from depth 2; to the surface of the 
disk. Hence 

^max 

/, = g S,{z) (1 - (9) 

Zrnin 

The model line profiles are calculated by convolving the 
disk emission with a gaussian beam, with the model 
beam size matched to the observations. 

4.1. Model Results - Axisymmetrie rotating disk 

A reasonable match to the observed profiles for the 
eastern nucleus is provided by a model in which the 
disk has center-to-edge line center optical depth ro,z = 4, 
dust optical depth ro,d = 1, velocity dispersion a = 85 
km s“^, and the rotation curve is flat with projected cir¬ 
cular velocity W = 135 km s“^. The normalizations of 
the line and continuum source functions are arbitrary; 
the continuum source function and the line center opti¬ 
cal depth have been adjusted to produce about the right 
levels for line and continuum for the disk center position 
and to match the absorption depth. In order to produce 
the deep central absorption the line center optical depth 
must be substantial (ro,z at least ^ a few) and there must 
be a substantial gradient in the CO J = 6 — 5 line source 
function (i.e., the excitation temperature) as a function 
of radius, decreasing outward. Thus, the deep absorption 
at the center is primarily from self-absorption. However, 
flat rotation curve models generically predict that the 
line peaks for the center position (p = 0) lie closer to 
the systemic velocity than those for larger offsets, in dis¬ 
agreement with the observations. A substantial improve¬ 
ment results from including a velocity perturbation (or 
a bump) in the rotation curve in the inner i? ^ 100 pc 
or so (compare the upper left panels of Figure [3 with 
the models shown in Appendix A). We have adopted the 
rotation curve shown in Figure El in which W peaks at 
about 240 km s“^. (The dip interior to the peak is just 
there to force V = 0 at R = 0; we have no sensitivity 
to emission (or its absence) on such small spatial scales.) 
The relative positions in velocity space of the red-wing 
peaks are in much better agreement with the observa¬ 
tions in the model which includes the velocity pertur¬ 
bation shown in the upper left panel of Figure 0 than 
for the flat rotation curve model shown in Appendix A 
(upper left panel of Figure [T5|). The velocity bump also 
has the effect of broadening the line profile somewhat, 
which also improves the match with the data0 We defer 
discussion of the observed rotation curve to §4.4. 

The exact shape of the curve should not be taken too 
seriously, as the needed V{R) depends on the spatial dis¬ 
tribution of the CO line emissivity, and non-circular ve¬ 
locities may also play a role. It is evident, however, that 
the velocity perturbation needs to be quite significant 
(^ 100 km s“^) — this is largely the result of the pro¬ 
jected velocity ^ 0 as the impact parameter p ^ 0. 
In the context of our models, it implies that there is a 
significant increase in the concentration of mass interior 
to R ^ 100 pc. 

® Although it is more subtle, this also improves the match be¬ 
tween the model and the data for the 0b/‘^ offset positions, for the 
lower-amplitude peaks. 






Eastern Nucleus - Major Axis 




Figure 5. The velocity profiles along the major axes of the Eastern (for PA ~35°) and Western (for PA '^111°) nuclei. The blue and 
magenta colored profiles are at the location of the line peaks and continuum peaks, respectively. The vertical line marks the systemic 
velocity of the system at 5350 km s“^. 



Radius (pc) 


Figure 6. Rotation curve adopted in our kinematic modeling of 
the eastern nucleus. Note: The dip interior to the peak is just 
there to force V = 0 at R = 0; our data do not have sensitivity to 
emission (or its absence) on such small spatial scales. 


The model and observed profiles are shown in Figure 
[71 The “red” and “blue” sides refer to which side of the 
double-sided profile dominates, so 1/ > 0 km s~^ (model) 
or "F > 5350 km s“^ (observed) for the red side — this 
is roughly the NE side of the disk (top panel), hence, 
the blue side is the SW side of the disk (bottom panel). 
The black line is for the central position (the continuum 
peak), while the yellow, orange, and red (green, cyan, 
blue) lines are for beam offsets ^b/2, Ob^ and 1.5^^ along 
the major axis for the red (blue) sides of the disk. 

The model makes a number of generic predictions that 
can be usefully compared with the observed spectra along 
the disk major axis: 

1. The disk center is where we see a symmetric, 
double-peaked line profile. 

2. The blue/red line peak brightnesses occur away 
from the center position. 

3. The blue/red line peak brightnesses and the extent 
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Figure 7. Comparison between the velocity profiles generated 
from an axisymmetric disk model and observations for the NE (top 
panels) and SW (bottom panels) sides of the eastern nucleus disk 
along the major axis of rotation. Black, yellow, orange and red 
(green, cyan, and blue) are for offsets of 0, 0^/2, Ob, and 1.5^^ 
away from the continuum center along the major axis approxi¬ 
mately towards the NE (SW) side of the disk. 

of the disk should be symmetric about the disk 
center. 

4. Emission in the blue wing is still present at a low 
level as one moves redward of the disk center, and 
vice versa. 

5. Producing line profiles that are as broad and ap¬ 
proximately Gaussian as observed requires a line- 
of-sight velocity dispersion that is comparable to 
the rotation velocity. See the appendix for a dis¬ 
cussion of what constraints the line profiles place 
on cr/ldot- 

Prediction (1) immediately implies that the disk center 
is at the continuum peak, not at the peak of the line 
emission, i.e., even though we find a 100 pc shift between 
the line peak and the continuum peak, the center of the 
disk is located at the continuum peak. This is what we 
would expect, as the dust emission is unaffected by the 
disk kinematics. It also implies that the central velocity 
of the deep absorption dip is the systemic velocity. 

Prediction (2), which is the inevitable result for a ro¬ 
tating disk with a source function gradient when the op¬ 
tical depth is not extremely high along all lines of sight, 
also agrees with the observations, but more so on the 
red (NE) side of the disk than the blue (SW) side. The 
disagreement becomes much more pronounced with the 
third prediction: the blue line peak maximum is consid¬ 
erably weaker than the red peak, and the disk extent 
on the blue side is noticeably less than that on the red 
side. EigureOshows that, while the ordering of peak flux 
density with position is basically correct, the observed 


line peak values are too small by about a factor of two 
compared to the model. In addition, the observed line 
widths of the blueshifted lines are systematically smaller 
than those of the redshifted lines, also by about a factor 
of two. 



Velocity (km/s) V^gdio (km/s) 



Velocity (km/s) (km/s) 


Figure 8. Comparison of the eastern nucleus kinematics between 
the axisymmetric disk plus foreground absorption model and ob¬ 
servations for the SW and NE sides of the nuclear disk along the 
major axis of rotation. Color coding as in Figure [7| 


The discrepancy with the axisymmetric model could be 
the result simply of asymmetry in the disk, in excitation, 
optical depth (e.g., the disk column density), or local 
kinematics — the disk dynamical timescale (radius of 
100 — 200 pc, rotation velocity ^ 130 km s“^) is a few 
million years. However, both the eastern and western 
nuclei show asymmetry in the same sense: although the 
CO emission is more extended than the continuum, the 
line emission in both nuclei is systematically weaker on 
the southern side of the disks (which is the blue side for 
the eastern nucleus), compared to the continuum extent 
(see Eigure[2fa,e,f)). 

There are two candidate mechanisms for producing 
such asymmetries in the line emission, absorption and 
outflow. We consider these in turn in the following sub¬ 
sections. 

4.2. Rotating disk models with foreground absorption 

The commonality of the CO line asymmetry in both 
nuclei suggests that something external to the nuclei is 
the cause. We therefore model it as additional, fore¬ 
ground absorption of the disk emission. This model as¬ 
sumes a Gaussian absorption profile in velocity charac¬ 
terized by the line center velocity Vabs, velocity disper¬ 
sion (Tabs^ and line center optical depth which are 
all allowed to vary with position. In practice, aabs = 70 
km s“^ provided a match in all cases. Because the 
















10 


line/continuum level is very important in determining 
the optical depth, we tweaked the model continuum lev¬ 
els slightly to precisely match the observations. 

Consider first the blue side of the disk, where we 
see large discrepancies between the original unabsorbed 
model and the data for blue velocities. The absorbed 
model emission profiles (after continuum subtraction) are 
shown in Figure [HI Given the simplicity of the model, the 
agreement is very good. Interestingly, there are signs of 
gradients as a function of position in both the absorption 
velocity and (clearly) in optical depth: the former shifts 
further away from zero velocity by about 40 km s“^, 
while the latter increases by nearly a factor of 3 (tq = 0.5, 
^- 0 b /2 = 0-9, T-Ob = 1-3, ^-30b/ 2 = 1-5). On the red 
side of the disk, only very weak absorption is allowed at 
P = with r no more than ^0.1, and basically none 

at larger offsets. (Indeed, at p = Ob^ there is clearly emis¬ 
sion at these velocities that is not present in the model.) 
Hence in this model there is little to no absorption on the 
red (NE) side, but the absorption optical depth is signif¬ 
icant at the continuum peak and it increases rapidly to 
the SW with distance from the center. 

Where might this absorbing gas he? One possibility is 
that both absorption features, i.e., the absorption that 
reduces the blue emission peak and the deep redshifted 
absorption (seen above 5505 km/s) are produced in one 
more-or-less contiguous feature. This would suggest that 
there is a very broad (^650 km s“^ in velocity space) 
complex of absorbing gas, the bulk of which is flowing 
outwards (blueshifted with respect to the systemic ve¬ 
locity), assuming that optical depth (which is probably 
larger for the blueshifted absorption than for the red- 
shifted) traces column density. The lower mass tail of 
absorbing gas extending to redshifted velocities may rep¬ 
resent dynamical spraying of the absorbing gas over a 
wide range of velocities, e.g., by tidal or ram-pressure 
stripping. 

To explain the common CO asymmetry and redshifted 
absorption of the nuclei, this gas must be at least as 
spatially extended as the nuclear separation (^400 pc). 
Given the complicated dynamics of the merging galaxies, 
the existence of such a spatial and kinematic feature is 
certainly plausible. Furthermore, the spatial location of 
the blue absorber would be coincident with the location 
of the redshifted absorption (Figures |2ff)), suggesting a 
physical connection between the two. 

Alternatively, the blue-shifted absorbing gas caus¬ 
ing the apparent line/continuum peak shift and the 
blue/red disk asymmetry could also be explained by 
an outer, kilop a rsec-s cale molecular disk , inferr ed by 
[Sakamoto et ah! (j2QQ9[ ) and [Scoville et al.[ (jl997[ ) from 
low-J CO interferometric observations. If this outer disk 
is suitably tilted with respect to the two embedded nu¬ 
clear disks, then the blue sides of the nuclear disks will 
appear fainter than the red sides as a result of absorp¬ 
tion. The magnitude of the required velocity dispersion 
(^ 70 km s“^) is certainly reasonable for such a disk, es¬ 
pecially as it has quite likely not settled into dynamical 
equilibrium. 

4.3. Rotating disk models with outflow 

[Sakamoto et al.[ (|2QQ9[ ) suggested that outflows with 
velocity Vout r\j 100 km s ^ are present in the nuclei of 


Arp 220, based on their reported detection of P Cygni 
profiles in the HCO+ (4-3) and (3-2) lines towards both 
nuclei and in the CO J = 3 — 2 line for the eastern 
nucleus, in SMA data with 0.3" resolution. Rotation 
is also clearly present in their data, and their preferred 
model included both rotation and outflow. 

With their more limited data set, they concluded that 
the absorption dip seen in the CO and HCO+ lines from 
the eastern nucleus is blueshifted by ^ 100 km s“^ with 
respect to the systemic velocity, which they estimated to 
be Vsys = 5415 ± 15 km s“^; a blueshift of the absorp¬ 
tion must be present if an outflow is responsible for the 
line profile asymmetries. In Figure [9] we overplot their 
line profiles on ours, for both the eastern (left panel) and 
western (right panel) nucleus. For the eastern nucleus it 
is clear that the velocities of the emission peaks and of 
the absorption dip for all three transitions are consistent, 
although the degree of asymmetry differs. We defer dis¬ 
cussion of the velocity shift of the absorption relative to 
the systemic velocity to §4.4. 

Since rotation is undoubtedly present, and the param¬ 
eter space for rotation plus outflow models is large, we 
adopted the same axisymmetric rotating disk model as 
in §4.1. As models in which there is a gradient in outflow 
velocity, rather than a constant velocity, provided better 
matches, we assumed a linear increase in outflow veloc¬ 
ity. (Models with exponential instead of linear rises differ 
insignificantly). The model outflow is symmetric about 
the line of sight and has a half-opening angle Oout- In 
practice, we found that models with Bout ^ 30° were in¬ 
distinguishable from models with larger half-angles, and 
provided better fits than narrower outflows. 

Reasonable approximations to the observed profiles are 
produced by models with Vout — 100 km s“^. The re¬ 
sulting proflles are compared with the observed profiles 
in Figure [TOl In detail the rotation plus outflow models 
do not match the data as well as the rotation plus ab¬ 
sorption models shown in Figure [HI especially on the blue 
side of the disk, where they tend to produce too much 
emission in the red wing of the lines. However, this con¬ 
clusion should not be considered in any way definitive, 
since the amount of unexplored parameter space is very 
large and we did not vary the underlying axisymmetric 
rotating disk model, only the outflow parameters. Hence 
it is only fair to conclude that the rotation plus outflow 
models for the eastern nucleus also provide an acceptable 
explanation for the observed line proflles. 

4.4. The Observed and Aetual Disk Veloeity Curves 

We extract the observed rotation curve for the eastern 
nucleus from the moment-1 map shown in Figure EJb) 
by taking a cut along the major axis. This is shown in 
Figure [TTl The curve has been color-coded as red- and 
blue-shifted with respect to the systemic velocity, which 
has here been defined to be the mean velocity at the 
position of the continuum peak, which is 5450 km s“^. 
To be specific, what is plotted here — the quantity cal¬ 
culated by moment-1 — is the intensity-weighted mean 
velocity (IWMV) for each pixel along the major axis. 
This weighting has an important influence on the rela¬ 
tion between the observed and actual rotation curves. 
Also, since the uncertainty in the mean velocity curve 
is completely dominated by the uncertainty in the posi¬ 
tion angle of the major axis, we have overplotted IWMV 
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Figure 9. Comparison of CO J = 6 — 5 line profile (this work) with the HCO+ (4-3)/(3-2) lines profiles from ISakamoto et al.l (|2009l i. 
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Figure 10. Comparison of the eastern nucleus kinematics be¬ 
tween the axisymmetric disk plus outflow model and observations 
for the SW and NE sides of the nuclear disk along the major axis 
of rotation. Color coding as in Figure [7] 

curves (with similar color-coding) for position angles be¬ 
tween 28° and 36°, which covers the range of uncertainty 
in the major axis position angle based on the continuum 
emission (which is unaffected by absorption) presented 
in Paper 1. The variation between them is insignificant. 

As mentioned in §4.1, the axisymmetric rotating disk 
model predicts that the continuum peak is the true cen¬ 
ter of the nucleus, and the central velocity of the self¬ 
absorption dip is the systemic velocity. This velocity is 
5350 km s“^ (e.g.. Figure [2fa)). Hence there is a dif¬ 
ference of approximately 100 km s“^ from the IWMV 
at the continuum peak. In addition, the observed rota¬ 
tion curve is noticeably asymmetric about the continuum 
peak IWMV. The red side is much shallower than the 
blue side: at 100 pc distance from the continuum peak, 
the velocity has risen by less than 100 km s“^ with re¬ 
spect to systemic velocity, whereas on the blue side it is 


over 200 km s“^. 

However, the asymmetry of the observed line profiles, 
and in particular the difference in profiles between the 
red and blue sides of the disk, produces a systematic off¬ 
set between the IWMV and the true velocity centroids 
of the emission. In the presence of such line asymme¬ 
tries (whatever their origin), IWMV curves can be very 
misleading about the underlying kinematics. As we now 
show, the intensity weighting of an asymmetric line pro¬ 
file can result in an IWMV curve which has both an offset 
from the true systemic velocity and an asymmetry in the 
derived rotation curve about the disk center, even though 
the underlying kinematic model is purely axisymmetric. 

Figure El (top) shows the IWMV curve produced by 
the axisymmetric model of §4.1, in the absence of ab¬ 
sorption (solid line). This is color-coded as in Figure 
nn and is of course symmetric about the disk center, 
since the line profiles, although asymmetric except for 
p = 0, are mirror images of one another on the red and 
blue sides of the disk. We define the apparent systemic 
velocity to be the IWMV at zero offset from the contin¬ 
uum peak (identical to the disk center), just as in Figure 
HB for the unabsorbed model this is of course equal to 
the true systemic velocity of 0 km/s. Q The dashed line 
shows the same thing, but now for the absorbed model of 
§4.2. On the red side of the disk, the absorption rapidly 
goes to zero with increasing distance from the disk cen¬ 
ter, so the red with-absorption curve converges with the 
no-absorption curve. Since the effect of absorption is to 
reduce the blue wing of emission, the absorbed IWMV 
curve shifts closer to zero on the blue side, and away 
from zero on the red side, with the result that the ab¬ 
sorbed IWMV curve shifts upward with respect to the 
unabsorbed IWMV curve. 

However, the inferred systemic velocity is no longer 
zero; it has been offset to the red, with Vsys — 125 
km s“^. (Recall that for the observed disk, the veloc¬ 
ity offset between the absorption dip central velocity, 
which should be the systemic velocity, and the IWMV 
at the continuum peak is AV ^ 100 km s“^ to the red.) 
In the lower panel we have subtracted off this apparent 

^ It is worth noting the weak resemblance between this “rotation 
curve” and the true disk rotation curve shown in Figure [7| This 
is a consequence of several systematic biases, especially the cos 6 
weighting of the line-of-sight velocity component, the relatively 
poor spatial resolution of the disk, the optical thickness of the 
lines, and the intensity-weighting used to derive the mean velocity. 
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Vsys from the absorbed model curve. This curve now ex¬ 
hibits all of the features shown by the observed IWMV 
curve: a steeper rise (and greater range) on the blue side, 
and a shallower red IWMV curve. Hence the observed 
IWMV curve is completely consistent with an underly¬ 
ing axisymmetric disk model, once the effects of the line 
profile asymmetries are taken into account. 



Figure 11. Observed IWMV rotation curve along the major axis 
of the eastern nucleus. The apparent systemic velocity is defined 
to be the IMWV at the continuum peak, 5450 km s“^. We have 
subtracted this quantity so that the plotted curves show rotation 
about the apparent systemic velocity, with the approaching and 
receding sides of the disk color-coded appropriately. Distances are 
positive in the NE direction (the red side of the disk). Rotation 
curves for position angles between 28° and 36° are overplotted to 
show that the variations between them are extremely small. 

In Figure [13] we show the same pair of IWMV curves 
as in Figure [12] for the outflow model of §4.3. Although 
the overall effect of the outflow is similar to that for ab¬ 
sorption, the apparent red side/blue side IWMV curves 
are noticeably less asymmetric than for the absorbed 
model, or as seen in the data. This is because the out¬ 
flow inevitably shifts line emission blueward from where 
it would occur in the absence of the outflow: as seen 
in Figure [TO] the absorption dip has shifted by about 
75 km s“^ to the blue from the true systemic velocity, 
and so part of the red wing of the line profile has also 
shifted to the blue, something which does not happen in 
the absorption model. In consequence, the IWMV at the 
continuum peak has shifted to the red by only about half 
as much as in the absorption model (about 68 km s“^ vs 
125 km s“^), and the two sides of the IWMV curves are 
more symmetric and less similar to the observed IMWV. 

While not a definitive argument, this does suggest that 
outflow models are less capable of matching the observed 
kinematic data than the absorption models. 

4.5. Rotation plus foreground absorption or rotation 
plus outflow? 

The line profiles and kinematics of the eastern disk are 
fit about equally well by a rotating, axisymmetric disk 
model combined with either foreground absorption on 
the blue side of the disk or with an outflow. On balance, 
however, we favor absorption as the explanation for the 
observed asymmetries in the line profiles and the disk 
extent with respect to the continuum peak: 

1. There is unequivocally foreground absorption 
present at redshifted (with respect to systemic) ve- 
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Figure 12. Rotation plus fore grou nd absorption: Model IWMV 
curves, color-coded as in Figure [TT] Top: Axisymmetric rotating 
model IWMV curve (solid line); with absorption (dashed line). 
The apparent systemic velocity is defined (as in Figure II) as the 
IMWV at zero offset from the disk center. With no absorption this 
is at the true systemic velocity, 0 km/s. The effect of absorption 
is to shift this velocity to the red by ~ 125 km s“^. Bottom: As 
the top panel, except that now the apparent systemic velocity has 
been subtracted from the dashed curve. Bottom: As the top panel, 
except that now the apparent systemic velocity (the IWMV at zero 
offset from the disk center) has been subtracted from the dashed 
curve. 

locities, over at least a few hundred km s“^, and it 
occurs spatially precisely where the blue shifted ab¬ 
sorption is expected to lie. Given the complicated 
dynamics and energetic nuclei present in Arp 220, 
the presence of additional absorption at blueshifted 
velocities (representing outflowing rather than in¬ 
falling gas) would not be unexpected. 

2. The line asymmetries produced by an outflow de¬ 
pend rather sensitively on the magnitude of the 
projected velocity compared to other velocities in 
the system, such as rotation. However, even though 
the two nuclear disks have very different appar¬ 
ent orientations, the asymmetries between the inte¬ 
grated line flux and the continuum emission for the 
two disks is very similar, as seen on the sky. This 
is natural in the absorption picture (the absorber 
only needs to be at least as large as the nuclear 
separation, and simply absorbs background radia¬ 
tion over its velocity extent) but rather contrived 
in the outflow model, where despite the different 
orientations of the nuclei, the projected rotation 
and outflow velocities conspire in such a way as to 
cause the line emission on the south sides of the 
disk to fade out with respect to the continuum in 
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Figure 13. Rota tion plus outflow: Model IWMV curves, color- 
coded as in Figure fm Top: Axisymmetric rotating model IWMV 
curve (solid line); with outflow (dashed line). Bottom: As the top 
panel, except that now the apparent systemic velocity (the IWMV 
at zero offset from the disk center) has been subtracted from the 
dashed curve. 

very similar fashion. (Recall that although we have 
not modeled the western nuclear disk, similar line 
profile asymmetries are seen there, e.g.. Figure [H) 

3. Despite the different orientations of the nuclei, the 
gas velocities on the south sides of the two nuclei 
(where the absorption occurs) are very similar (Fig- 
urelSfb)): from the eastern to the western nucleus 
on the blue side, all of the emission falls within 
the same range (covering ^100 km s“^) so the ab¬ 
sorbing gas can occupy the same range of velocities 
over its whole extent, with no need to introduce a 
velocity gradient across the absorber. 

4. The IWMV curve produced by the rotation plus 
absorption model resembles the observed curve 
much more closely than the rotation plus outflow 
model, in considerable part because in the outflow 
model some of the red side emission must be shifted 
blueward. 

The rotation plus absorption model still requires that 
there is outflowing, blueshifted gas in the Arp 220 sys¬ 
tem. However, this gas is not dynamically responsible 
for the line asymmetries. Whichever model is correct, 
however, the mass of gas involved in the outflow is very 
substantial. In the rotation plus outflow model, a sub¬ 
stantial fraction of the total mass of the disk (for which 
we estimate a warm gas mass 1.3 x IO^Mq, see 

Table 1) must be participating in the outflow in order to 


produce the observed line profile asymmetries, the exact 
amount depending on the opening angle of the outflow. 
This also implies that the lifetime of the disk is only 
tout ~ R/Vout ~ 10® yr. 

For the rotation plus absorber model, the mass is more 
uncertain. We can write the line center optical depth in 
a Doppler-broadened line as 


To = 
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where Ui is the line frequency, Aji is the Einstein A- 
coeflicient, the Qi and gj are the statistical weights of the 
levels, Ni is the column density in the lower state of the 
transition, and AV is the Doppler width. For the CO 
J = 6 — 5 line, Aji = 2.126 x 10“^ s“^ and n = 691.473 
GHz. For the blue-side absorber, AV = ^/2 x 70 km s“^ 
and To ^ 1. The H 2 column density corresponding to 
equation [To] is then 


( 11 ) 

where xqo = 10 ^ is the CO abundance with respect to 
H 2 and Qr = is the rotational partition 

function for CO. 

Spatially, the absorber must extend at least ^ 400 pc 
(the separation of the nuclei) by 220 pc (the extent of 
the observed absorption on the SW side of the eastern 
nuclear disk), giving a minimum area Aabs ~ 8-8 x 10^ 
pc^. The minimum mass of the absorber is then 


Mabs > 

?a4.2xl0®- - V .A l - ^ M(^^) 

(a;co/10-4) V J 

where /a is the mean molecular weight per H 2 molecule. 

The partition function Qr is almost certainly > 10, 
and plausibly as large as 100, while the correction for 
stimulated emission is likely ^ 1 to a few, giving 

Mabs ~ few X ( 10 ^ - 10®) Mq. (13) 

This is again a very substantial amount of mass, com¬ 
parable to the total mass of warm molecular gas in Arp 
220 (Note: The warm gas mass in Arp 220 is ^10% of 
the total gas mass). The velocity characterizing the ab¬ 
sorber is much larger than in the outflow model, reaching 
blueshifted velocities of nearly 300 km s“^ with respect 
to systemic, as this gas must absorb the blue wings of 
the emission lines. 

An interesting question is whether we should expect 
to see the blue absorber emission away from the nuclear 
disks. Exploratory non -LTE modeling using RADEX 
(|van der Tak et ahlDOOTl) of plausible absorber parame¬ 
ters (defined as those which produce Qr not much larger 
than 10, so as to minimize the mass of the absorber) 
suggests that the ratio of warm gas to absorber peak 
antenna temperatures in the J = 6 — 5 line is close to 
a factor of 10. Eurthermore the velocity extent of the 
absorber is almost certainly substantially smaller than 
that of the warm gas emission, further increasing the ra¬ 
tio in integrated intensity. Comparison with Eigure|2fa) 
suggests that it is unlikely that the absorber emission is 
detectable on top of the general nuclear disk emission. 














14 


There is clearly extended emission present in the chan¬ 
nel maps between 5105 km s“^ and 5225 km s“^, where 
we would expect the blue absorber emission to exist, but 
its impossible to conclude if this is from the absorbing 
gas rather than the general disk emission. 

ALMA observations of other CO transitions (e.g., J = 
4 — 3, 7 — 6) will be crucial in constraining the parti¬ 
tion function, and excitation temperature, by 
providing estimates of the variation of line optical depth 
with J, and thus reducing the uncertainty on the mass es¬ 
timate of the absorber. For example, the same absorber 
models described above predict that the absorber is con¬ 
siderably brighter in the J = 3 — 2 line, both absolutely 
and relative to the warm gas (cf. end of §3 and Figure 

m). 


5. DYNAMICAL MASS AND DISK STABILITY 


The numbers derived from the kinematic modeling in 
§4 can be used to estimate the enclosed mass in the east¬ 
ern nucleus. Adding the asymptotic rotation velocity 
(135 km s“^) and velocity dispersion (85 km s“^) in 
quadrature, we derive a dynamical mass within r = 220 
pc of Mtot ^ 1.6 — 2.9 X IO^Mq, where the larger number 
is obtained by correcting the rotation velocity (but not 
the velocity dispersion) f o r an in clination of 50°, obtained 
by iBarcos-Muhoz et al ] (IMl from a two-dimensional 
nonlinear least-squares fit of a thin, tilted, exponential 
disk to their 33 GHz radio continuum image with 0.08" 
X 0.06" resolution. Given the uncertainties, this enclosed 
mass is in agreement with the total molecular gas mass of 
1.1 X 10^ derived in Paper I from the continuum emission 
(adopting a gas-to-dust mass ratio of 100). This implies 
that the dynamical mass within the radius of 220 pc is 
largely dominated by the molecular gas, consistent with 
previous studies of Arp 220. 

An interesting question is whether the nuclear disks are 
liable to be gravitationally unstable to fragmentation and 
the formation of stars given the large inferred mass sur¬ 
face densities (see Paper I) of about 5.4 x 10^ Mq/pc^ and 
1.4 X 10^ Mq/pc^ associated with the disks of the east¬ 
ern and western nuclei, respectively. We investigate this 
using our model for the eastern n ucleus and the Toomre 
Q-parameter for a gaseous disk (jToomrd 1 19641) . in the 
form 




VsK. 

ttGS 


> 1 


(14) 


for stability, where Vg is the sound speed in the gas, k is 
the epicyclic frequency, G is the gravitational constant, 
and S is the gas surface density^ 

To evaluate S, we use our dynamical mass estimate 
Md ^ 2 X 10^ Mq and an effective disk area A = 7r0‘^, 
where Oe = 276 pc is the deconvolved FWHM for 
the eastern nucleus. This gives S ~ 8400 Mq pc“^. 
Note that this is an average surface density (and thus 
lower) compared to the value reported in paper I, which 
was effectively the peak surface density calculated on 
a smaller spatial scale (radius < 50 pc). We also as¬ 
sume a flat rotation curve with W = 135 km s“^, so that 
K, = ^/2Q = which we evaluate at = ^£;/2 as 


® Strictly speaking, this equation applies to a thin disk, which is 
not the case here, since the ratio of the velocity dispersion to the 
rotation velocity cr/Vc ~ 0.5 — 1. The effect of t his will be to make 
the disk more stable than implied by equation (d. 


a typical value. Finally, we use the inferred velocity dis¬ 
persion cr = 85 km s~^ as an effective sound speed. With 
these parameters, we find Q ^ 1.5; this number should be 
regarded as uncertain by at least a factor of two (further¬ 
more, S, H, and a are all functions of radius). However, 
the estimated Toomre Q-parameter for the gas disk of the 
eastern nucleus is of order unity, suggesting that the disk 
is marginally stabilized against gravitational fragmenta¬ 
tion and collapse by the large velocity dispersion. This 
large velocity dispersion requires energy input to sustain 
it unless we are catching Arp 220 in a very short-lived 
phase. The amount of random (non-rotational) kinetic 
energy in the ISM can be written 

Ekin ~ - 10 “ 

where M is the mass of molecular gas and aioo is the 
velocity dispersion in units of l QQkms~b Dim ensionally, 
the energy dissipation rate is (|Malonevlll999[) 

Ekin a 

~6 X ( io^q ) -^ioVioo erg 8“^ (16) 

where R = IQQi^ inn pc is a chara cteristic radius and the 
coefficient r] <1 ()Mac Lowlll999[ ). The ratio of these two 
equations gives a dissipation time 

tdiss = --5 X 10"^-yr. (17) 

Ekin Wioo 

(A very similar estimate is obtained simply from assum¬ 
ing that the areal filling factor is roughly unity and that 
collisions are strongly dissipative.) For the eastern nu¬ 
cleus of Arp 220, Rioo ^ cnoo ^ 1, which gives the dissi¬ 
pation time of only tdiss ^ 10^ yr. This short timescale 
implies that there must be an energy input sustaining 
the large velocity dispersion. 

6. WARM MOLECULAR GAS MASS AND LUMINOSITY 

We determined the mass of the warm molecular gas 
(listed in Table 1) in the individual nuclei using the in¬ 
formation f rom our Herschel SP IRE-FTS observations 
of Arp 220 (|Rangwala et al.|[2011[) . We derived the ratio 
of the total warm molecular gas mass to the total GO 
J = 6-5 line luminosity in Arp 220; the former is de¬ 
rived from non-LTE modeling of the full GO SLED from 
J = 1 — Oto J = 13 — 12 and the latter is measured 
directly from observations. This ratio M(H2)/I/co(6-5) 
is about 21 Mq/Lq. Using this ratio and the individ¬ 
ual integrated GO J = 6 — 5 ALMA fluxes for the two 
nuclei we estimate the warm molecular gas mass to be 
1.4 X 10^ Mq and 1.6 x 10^ Mq, in the eastern and west¬ 
ern nucleus, respect ively0 These are roughly 1/lOth of 
the total molecular gas masses measured from dust con¬ 
tinuum emission in Paper I and are consistent with the 
war m molecular gas mass fra ction found in nearby galax¬ 
ies (|Kamenetzkv et al']l2Q14D . 

® We note that the beam size of Herschel is large compared to 
the angular extent of Arp 220 and therefore this ratio is an aver¬ 
age quantity over the whole system and is likely to vary spatially. 
Hence, our mass estimate for the warm molecular gas for the indi¬ 
vidual nuclei should be used with caution. 
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We list the observed I/coj= 6 - 5 /^fir ratio in the two 
nuclei in Table 1 using the Tfir value calculated in Pa¬ 
per 1. This ratio differs by almost an order of magnitude 
between the two nuclei - a direct result of their Lfir val¬ 
ues differing by an order of magnitude. The ratio in the 
eastern nucleus of ^ 6.7 x 10“^ is much higher than the 
western nucleus; the value for the latter is within a factor 
of two o f the average value of ^ 10“^ found in nearby 
galaxies (jKamenetzkv et ^l2Q14f ) . Despite the disparity 
in the ratio of CO J = 6 — 5 luminosity to bolomet- 
ric luminosity, the line-to-continuum ratios at 691 GHz 
(listed in Table 1) are the same for the two nuclei, im¬ 
plying that the eastern nucleus contains an additional 
power source which is depositing energy into the gas 
without affecting the dust continuum emission. Such an 
energy source is consistent with mechanical energy from 
shocks. A similar scenario is inferred in another nearby 
IR-luminous galaxy, NGC 6240, in which the GO is also 
unusually luminous relative to the bolo metric luminos¬ 
ity, w ith Tcoj= 6 - 5 /^fir ^ 7.6 x 10“^ (jMeiierink et al.l 
l2Q13f ) . It is intriguing that even though there is evidence 
of large-scale gas dynamics (both outffow and inffow) to¬ 
wards both nuclei, only the eastern nucleus - which is 
the less luminous of the two, and which has not been 
suggested to host an AGN - exhibits enhanced CO lumi¬ 
nosity relative to its FIR luminosity. 

The large velocity dispersion found for the east¬ 
ern nucleus is one possi ble source of energy f or the 
warm CO emission. Fro m iRangwala et aP (j2Ql If ) (also, 
iKamenetzkv et al']l2Q14f) . the total CO luminosity of Arp 
220 is Leo ~ 2 X lO^L©. We assume that half of this 
luminosity comes from the eastern nucleus. However, at 
the warm temperature inferred for this component, most 
of t he gas cooling will b e from molecular hydrogen, no t 
CO (jNeufeld et al.llTo^ . From iLe' Hour lot et al.l (jl999f ). 
the H 2 cooling function is Ahs ~ 10“^^ erg s~^ (=11 
Lq/Mq) with a variation of about a factor of three over 
the range of temperatures and densities that characterize 
the warm molecular gas in Arp 220. With a warm gas 
mass Mh2 — 1.3 x 10^ Mq for the eastern nucleus (Table 
1), the total warm gas cooling Twarm ^ 6.5 x 10^^ erg 
s“^ (with a range of a factor of 3). Comparison with 
equation M shows the energy dissipation rate is of the 
right order of magnitude to excite the warm molecular 
emission. This new result disagrees with the conclusions 
of Rangwala et al (2011). Their much smaller estimate 
of the energy available from turbulent heating was based 
on a gas velocity gradient derived from modeling of the 
global CO emission, while our new value uses the directly 
observed velocity dispersion provided by our ALMA ob¬ 
servations. 

We also consider the energy available from super¬ 
novae and stellar winds. Using th e observed supernovae 
rate of z/sn = 4 ± 2 per year (jLonsdale et al.l l2006[) 
and (F'sn = 10 ^^ erg) in the following expression from 
iMaloneH (j 19991 ): 


r o in43/ ^SN x/ ^SN x _l 

Lsn^3x10 (-—rT)(T7;5i- )ergs , 

lyr ^ 10^^ ergs 


(18) 


we obtain a total mechanical energy of 10^4 g 
The energ y output from stellar wm ds is similar to su¬ 
pernovae (|McCrav fc Kafat^ 119871 ) , giving a total of 
^ 2 X 10^^ erg s“^, also large enough to power the 


warm gas emission and maintain the observed veloc- 
ity dispersion; this is consistent with the conclusions of 
IRangwala et al.l (j2QllF 

7. HIGHLY EXCITED SIO ABSORPTION TOWARDS THE 
WESTERN NUCLEUS 

Absorption from the SiO J = 16 — 15 line at 694.293 
GHz is seen towards the location of the continuum emis¬ 
sion from the western nucleus (Figure [2fe) and Figure [3]). 
The absorption is point-like with a full velocity width of 
at least 400 km s“^; the line width could be as large as 
500 km s“^ if the SiO absorption continues through the 
small gap in our spectral coverage. The mean depth of 
the absorption in the central 200 km s“^ of the line is 
-0.23 Jy beam“^ or -6.9 K. No significant SiO absorp¬ 
tion is seen towards the eastern nucleus; any absorption 
in this nucleus must be at least 12 times fainter over the 
same relative velocity range. 

The SiO J = 5 — 4 line at 217 GHz has been prev i- 
ously seen in Arp 220 in emission (jMartm et al.ll20lTl ). 
although the angular resolution of the dat a was insuffi¬ 
cient t o distinguish between the two nuclei. iMartm et al.l 
(I 20 T 1 I) note that the SiO emission is not particularly en- 
hance d in Arp 220 compyed to other species. In con¬ 
trast, [GarHiiBurniQ^ES] (j2010[ ) find that SiO emission 
near the AGN in NGG 1068 has an enhanced abundance 
by a factor of 10-100 con ipared to locations in the star- 
burst ring. The fact that iMartm et al.l (j2011f ) do not see 
an enhanced SiO abundance may be related to the in¬ 
ability of their data to isolate the region nearest to the 
AGN from surrounding starburst re gions. _ 

For Galactic star forming regions. ISchilke et al.l (jl997[ ) 
showed that SiO is released via neutral particles collid¬ 
ing with charged grains. Shock speeds of 25 km s“^ and 
densities of 1 x 10^ cm“^ are sufficient to reproduce the 
observed SiO column densities. The mean density of the 
dusty disk in the western nucleus of Arp 220 is estimated 
to be > 10^ cm“^ (Paper I) and high turbulent veloci¬ 
ties would be consistent with the large CO line-widths 
seen towards this nucleus. Thus, we might expect sig¬ 
nificant SiO to be present throughout this dense dusty 
disk. However, the SiO J=16-15 absorption is clearly 
more compact than the dusty disk; the strength of the 
SiO absorption does not change as the beam is increased 
to a FWHM of 0.5", while the peak continuum fiux in¬ 
creases from 1.15 to 1.69 Jy beam“^. Thus, high exci¬ 
tation SiO molecules appear to be concentrated towards 
the center of the western nucleus. However, because the 
dust emission is quite optically thick (r ^ 5, Wilson et 
al. 2014), the radius of the region producing the SiO ab¬ 
sorption has to be comparable to the scale height of the 
disk. 

The question arises as to how SiO is excited to such a 
high energy level in Arp 220. The upper state of the SiO 
J = 16 — 15 line is 283 K above ground and the transition 
has an Einstein coefficient logA^j = —1.74 (taken from 
the Splatalogue molecular database). Thus, the critical 
density of this transition is ^ 2 x 10^ cm“^, 1000 times 
higher than that of the CO J = 6 — 5 line. Because the 
SiO J = 16 — 15 line is seen in absorption, the excitation 
temperature of this transition must be < 100 — 200 K, 
the estimated brightness temperature of the continuum 
source at this frequency. One possible mechanism to ex¬ 
cite SiO is via radiative pumping. Radiative pumping 
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has been invoked to explai n the hi^h-J absorptio n lines 
of HCN seen in Arp 220 (jRangwala et alJl2Qllf ). The 
first excited vibrational state of SiO is at ^ 8 /im, near 
the peak of a 360 K blackbody spectrum. Given the 
high mean temperature of the western nucleus, radiative 
pumping by the intense radiation field of an AGN or an 
extreme nuclear starburst appears to be a viable mecha¬ 
nism. 

The alternative would be collisional excitation in the 
inner, presumably higher density regions, of the nuclear 
disk. However, a 50 pc radius disk with a mean density 
of 4 X 10^ cm ^ requires a steep density profile r to 
achieve densities comparable to the critical density before 
the dust reaches its sublimation temperature. Therefore, 
it is difficult to see how a sufficiently extended region of 
excited SiO could be achieved under these conditions. 

8. CONCLUSIONS 

Our Gycle-0 ALMA observations of the CO J = 6 — 5 
line show complex features in the morphology and kine¬ 
matics of the warm molecular gas not seen before, neither 
in low-J lines tracing cold gas nor in lower-resolution CO 
J = 6 — 5 observations from the SMA. The overall mor¬ 
phology of the warm emission is similar to that traced by 
the low-J lines. However, we detect an apparent offset 
of ^100 pc between the peak of the continuum emission 
and the peak of the line emission. This offset is present 
in both nuclei and is in the same direction. Compari¬ 
son with kinematic models suggests that the offset arises 
due to the absorption of the entire South (blue) side of 
the central region containing the two nuclei by gas ly¬ 
ing in front of the nuclei. There is also clear evidence 
in our data for redshifted absorption, which could be 
interpreted as an inflowing molecular filament. If a sin¬ 
gle coherent feature is responsible for all of the absorp¬ 
tion features - both the redshifted absorption and the 
blueshifted absorption that produces the apparent offset 
and the line profile asymmetries - it has a total veloc¬ 
ity extent of ^700 km s“^ and a spatial extent of ^400 
pc. Alternatively, the blueshifted absorption may arise 
in a larger (kpc)-scale molecular disk, as suggested by 
earlier interferometric observations, and have no dynam¬ 
ical connection with the infalling gas. Although similar 
asymmetries in the line profiles can be produced by a 
combination of rotation plus outflow rather than rota¬ 
tion plus absorption, on balance we believe that the data 
favor absorption rather than outflow as the source of the 
asymmetries. 

From our modeling, we conclude that the dynamical 
center of the eastern nuclear disk is at the location of 
the continuum peak, and that a large gradient in the 
excitation is needed to explain the deep absorption de¬ 
tected at the center of the eastern nucleus. In addition, 
the velocity dispersion in the disk must be a significant 
fraction of the rotation velocity to be consistent with the 
large line widths observed in the two nuclei. The dy¬ 
namical mass estimated from our kinematic model im¬ 
plies a very large mass surface density for the eastern 
nucleus, suggesting the disk could be gravitationally un¬ 
stable and prone to disruption. However, the estimated 
Toomre Q-parameter for the gas disk of the eastern nu¬ 
cleus is (marginally) greater than unity, suggesting that 
the disk is plausibly stabilized against gravitational frag¬ 
mentation and collapse by the large velocity dispersion. 


We measure an unusually high ratio of CO J = 6 — 5 
luminosity to the total FIR luminosity in the eastern 
nucleus compared to the western nucleus; for the latter 
the value is consistent with the average found in nearby 
galaxies. This suggests that there is an additional energy 
source, such as mechanical energy from shocks, present 
in the eastern nucleus that is generating more CO lumi¬ 
nosity without affecting the dust continuum emission and 
hence the total FIR luminosity. We find that the total 
warm gas cooling from H 2 approximately matches the 
energy dissipation rate (from velocity dispersion), sug¬ 
gesting that this energy dissipation could also potentially 
power the warm molecular emission, in addition to the 
energy from super novae and stellar winds. 

Without the exceptional sensitivity and spa¬ 
tial/spectral resolution from ALMA it would have 
not been possible to detect the complexity in the 
morphology and kinematics of the warm molecular 
gas in Arp 220, and to measure the various quantities 
associated with them with high signal-to-noise ratios. 
This is especially important for detecting the signature 
of inflowing molecular gas, which is much harder to 
detect. The high S/N of finely sampled velocity profiles 
allowed us to distinguish between various kinematic 
models. It will be highly beneficial to have the same 
quality of observations for the multiple CO lines that 
are accessible with ALMA, to look for spatial variations 
in the gas excitation and better characterize the mass 
and excitation of the large-scale gas dynamics (both 
outflow and inflow) towards the two nuclei. This kind 
of spatial/kinematic analysis is highly complementary 
to the analysis of Herschel observations described in 
§1. Our accepted Cycle-2 observations of Arp 220 (and 
NGC 6240) will provide a step forward in this direction. 
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17 


APPENDIX 

APPENDIX A: ADDITIONAL MODELING RESULTS 




Velocity (km/s) Velocity (km/s) 

Figure 14. Velocity profiles generated from an axisymmetric disk model for two cases: (Top) Doppler parameter (\/2 x dispersion) = 
180 km s“^ and Vmax = 0 km s“^ and (Bottom) Doppler Parameter = 10 km s“^ and Vmax = 180 km s“^. Color coding as in Figure [7| 






Figure 15. Model profiles for different values of Doppler parameter and maximum rotation velocity with their quadrature sum fixed. 
Color coding as in Figure [7| 

The strength of the emission on the opposite-velocity side of the absorption dip as one moves away from the center 
along the major axis towards either the red or blue sides of the disk, the relative peak heights along the major axis, and 
the shapes and velocity extents of the red and blue line profiles, can be used as diagnostics of the relative magnitudes 
of the velocity dispersion and rotation velocity. Since there are clearly systematic issues with the blue side of the 
disk, discussed in Section 4, we mostly consider the constraints provided by the red side of the emission. In the limit 
of zero rotation velocity (assuming the same optical depth and source function parameters as in Figure [7j), the line 
profile is symmetric about the absorption dip, and the line peak flux densities are greatest at the disk center, while if 
the velocity dispersion is much smaller than the rotation velocity, the line profiles are very non-Gaussian, and peak 
sharply at the maximum rotation velocity at all positions along the major axis (see Figure d. To illustrate the 
approximate allowed ranges of the rotation velocity and velocity dispersion, in Figure [15] we plot the resulting line 
profiles (at half-beamwidth intervals and Gaussian-weighted, as in Figures [7]) for four different models, each labeled by 
Vc and AVd- In all cases Vc and AVn have been varied so that their quadrature sum is nearly fixed, which ensures that 
the peak-to-peak separation of the red and blue peaks is similar to the observed (~ 350 km s“^) in all models. Note 
that this separation is generally significantly larger than twice the rotation velocity, as is the extent of the emission 
in velocity space: this is a consequence of both the substantial velocity dispersion and of optical depth effects, as the 
models are not optically thin. For clarity, we only plot the models for the center and red-side positions; the blue-side 
emission lines are the mirror-image of the red-side. 
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There is little change in the peak flux densities as AVd is increased to 135 km s“^; the main change is a modest 
increase in the line widths. As is evident in Figure [3 the model line profiles are a little narrower than the observed, 
and do not exhibit the trend of increasing linewidth with distance along the major axis seen at the Ob and 1.50b offset 
points, indicating deviations from the model assumptions (constant Vc and AVd)- It does predict a slightly stronger 
red peak at AO = —ObI‘^ (the blue side of the major axis) than the AVd = 120 km s“^ model, which is already slightly 
too strong compared to the observations. 

Raising the rotation velocity to 155 km s“^ makes the lines narrower, slightly reduces the difference in peak flux 
density between the center and AO = 30bI positions, and substantially reduces the intensity of the blue-shifted emis¬ 
sion. Reducing Vc to 90 km s“^ and raising AVd correspondingly (i.e., swapping the values from the previous model) 
makes the center and AO = Ob peak flux densities nearly equal, in substantial disagreement with the observations, 
and makes the blue-shifted emission substantially brighter. Finally, raising Vc as high as 175 km s“^ produces nearly 
equal peak flux densities in the center and AO = 36>b/ 2 positions, and more nearly equal peak flux densities in the 
AO = Ob 1^^ and AO = Ob I positions, again in disagreement with the observations. 

In summary, these model variations suggest that the disk rotation velocity is unlikely to be substantially below 120 
km s“^, or above 155 km s“^, with the velocity dispersion scaled accordingly. 
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